基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

涡轮基组合循环发动机分布式控制系统通信网络混合拓扑结构优化方法

涡轮机组合循环(Turbine based combined cycle,TBCC)发动机控制系统通信网络拓扑结构是其分布式控制系统方案设计的重要部分,优化网络拓扑结构可提高发动机推重比和控制系统可靠性。本文基于智能优化算法提出TBCC分布式控制系统网络拓扑结构优化方法。基于图论建立TBCC几何模型和网格模型,以重量和可靠性为优化性能指标,同时考虑发动机表面高温区域以及控制节点的工作可靠性,分别采用粒子群算法和遗传算法优化星形结构中智能中央节点位置、中央节点的环形拓扑结构,获得星形-环形混合拓扑结构。仿真实例表明,基于本文方法优化所得的混合拓扑结构相较于星形集中式控制结构,系统重量降低了51.9%。     

面向多属性条件的空间对象检索算法

针对现有空间索引不能满足多样化的检索需求,提出两种新型空间索引,能够同时面向空间属性、文本属性与划分属性进行空间对象检索。将分类技术应用于空间对象检索中,提出了基于划分索引与IR-Tree的混合索引以及先划分索引再IR-Tree的索引方法,不仅满足了多样化的空间检索需求,而且有效地解决了传统空间索引更新维护代价大的问题。基于真实的北京市POI数据集进行实验,结果表明本文提出两种索引是有效的且检索效率高。与传统空间索引相比,提出的空间索引很好地解决了具有划分属性的空间对象检索问题,并且具有较高的检索效率。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

一种地理围栏内无人机航迹跟踪方法

针对UTM体制中无人机在地理围栏内的飞行监视问题,提出一种约束状态相关模态转换混合估计算法(CSDTHE)。采用随机线性混杂系统模型对无人机运动状态进行建模,利用CV、CT和CA三种模态描述无人机的飞行状态,以构建地理围栏内无人机运行的通用模态转换模型框架。利用飞行模态改变点(FMCP)定义相关模态转换参数,设计模态转换条件,生成模态转换概率矩阵,从而建立与状态相关的模态转换模型。运用约束卡尔曼滤波(CKF)方法对直线阶段和转弯阶段的无人机运动速度分别施加等式约束,并通过仿真实验验证了CSDTHE算法对无人机跟踪的有效性。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

不同迎角下脊形前体绕流数值模拟研究

脊形前体有较强的背风涡流场,不同的前体形状对前体涡流场和气动力有很大的影响。本文针对脊形前体飞行器大迎角湍流大分离流动计算的困难,采用IDDES混合湍流模型,以及与之匹配的非定常算法,研究了不同来流迎角下脊形前体的气动特性,以及背风涡非定常演化、破裂的细致流动结构。选取了不同脊形角,以及不同上、下高宽比的脊形前体进行计算。计算结果表明,在迎角较小时,随着迎角的增大,前体主涡会逐渐增强,在迎角较大时,前体主涡破裂;在相同迎角下,脊形角较小时,前体涡较强,涡升力也更大;对于相同脊形角的前体,当上半截面高宽比较小时,前体主涡强度较大,前体涡破裂临界迎角较小,即会提前破裂。     

Multi-scale design and optimization for solid-lattice hybrid structures and their application to aerospace vehicle components

By integrating topology optimization and lattice-based optimization, a novel multi-scale design method is proposed to create solid-lattice hybrid structures and thus to improve the mechanical performance as well as reduce the structural weight. To achieve this purpose, a two-step procedure is developed to design and optimize the innovative structures. Initially, the classical topology optimization is utilized to find the optimal material layout and primary load carrying paths. Afterwards, the solid-lattice hybrid structures are reconstructed using the finite element mesh based modeling method. And lattice-based optimization is performed to obtain the optimal cross-section area of the lattice structures. Finally, two typical aerospace structures are optimized to demonstrate the effectiveness of the proposed optimization framework. The numerical results are quite encouraging since the solid-lattice hybrid structures obtained by the presented approach show remarkably improved performance when compared with traditional designs.     

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

云制造环境下机床装备资源选择方法

针对云制造环境下机床装备资源具有数量大、异质异构、服务制约因素多等问题，提出了基于多准则决策的机床装备资源选择方法。首先，利用本体论对机床装备资源建模，将物理资源虚拟化封装，实现虚拟资源与物理资源的逻辑映射。进一步地，使用语言变量量化决策者评估信息，采用模糊决策与试验评价实验室（FDEMATEL）和熵值法（EW）相结合的组合赋权法，引入模糊VIKOR（FVIKOR）多属性决策方法对机床装备资源进行综合评估，建立FDEMATEL-EW-FVIKOR的多准则混合决策模型。最后，将该方法应用于某生产制造企业，实验对比决策结果验证了多准则混合决策模型的可行性和有效性。